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AERONAUTIC  SYMBOLS 
1.  FUNDAMENTAL  AND  DEKIVED  UNITS 


Symbol 

Metric 

English 

Unit 

Abbrevia¬ 

tion 

Unit 

Abbrevia¬ 

tion 

Length _ — 

Time - 

Force . 

l 

t 

F 

meter _ 

second _ 

weight  of  1  kilogram - 

m 

8 

kg 

foot  (or  mile) - 

second  (or  hour) - 

weight  of  1  pound - 

ft  (or  mi) 
sec  (or  hr) 

lb 

Power _ — 

Speed - 

P 

V 

horsepower  (metric) - 

r kilometers  per  hour . 

\meters  per  second - 

kph 
mps  ... 

i 

horsenowor  . . — 

1  miles  pcr  hour . -- 

feet  per  second - 

h‘\ 

mph 

fps 

2.  GENERAL  SYMBOLS 


Weight  =  mg 

Standard  acceleration  of  gravity=9.80665  m/s 
or  32.1740  ft/sec* 

W 

U&ss=j 

Moment  of  inertia=mP.  (Indicate  axis  of 
radius  of  gyration  k  by  proper  subscript.) 

Coefficient  of  viscosity 

,  8.  AERODYNAMIC  SYMBOLS 


v  Kinematic  viscosity 

n  Density  (mass  per  unit  volume)  _ 

Standard  density  of  dry  air,  0.12497  kg-m  -s  at  lo  C 

and  760  mm;  or  0.002378  lb-ft‘4  sec3 

Specific  weight  of  “standard”  air,  1.2255  kg/m  or 
0.07651  lb/cu  ft 


Area 

Area  of  wing 
Gap 
Span 
Chord 

Aspect  ratio,  g 
True  air  speed 
Dynamic  pressure,  ^pT73 

Lift,  absolute  coefficient  CL— 

Drag,  absolute  coefficient  CD 


L 

qS 

“2  S 


Profile  drag,  absolute  coefficient  0D 
Induced  drag,  absolute  coefficient  CD 


D0 

qS 

J> t 
'*  qS 

Dv 


Parasite  drag,  absolute  coefficient  CDp— 
Cross-wind  force,  absolute  coefficient  C  c '  “g 


a 


Q 

c 

E 


a 

e 

a0 

«< 

aa 


Angle  of  setting  of  wings  (relative  to  thrust  line) 
Angie  of  stabilizer  s riling  (relative  to  thrust 

line) 

Resultant  moment 
Resultant  angular  velocity 

Reynolds  number,  p-J-  where  l  is  a  linear  dimen¬ 
sion  (e.g.,  for  an  airfoil  of  1 .0  ft  chord,  100  mph, 
standard  pressure  at  15°  C,  the  corresponding 
Reynolds  number  is  935,100;  or  for  an  airfoil 
of  1.0  m  chord,  100  mps,  the  corresponding 
Reynolds  number  is  6,865,000). 

Angle  of  attack 

Angle  of  downwash 

Angle  of  attack,  infinite  aspect  ratio 

Angle  of  attack,  induced 

Angle  of  attack,  absolute  (measured  from  zero- 
lift  position) 

Fliglit-path  angle 
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By  Dean  R.  Chapman 

Page  2,  equation  (6),  the  exponent  32  should  read  3/2. 

Page  2 9  line  preceding  equation  (7),  "than"  should  read  "then, 
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LAMINAR  MIXING  OF  A  COMPRESSIBLE  FLUID 

By  Dean  It.  Chapman 


SUMMARY 

A  theoretical  invest ujat ion  of  the  velocity  profile*  jor  (am mar 
mixing  of  a  high-velocity  stream  with  a  region  of  fluid  at  rest 
has  been  made  assuming  that  tht  Prandtl  number  is  unity. 
.1  method  which  involves  only  quadratures  is  presented  for 
calculating  the  velocity  profile  in  the  mixing  layer  for  an  arbitrary 
value  of  the  free-stream  Mach  number. 

Detailed  velocity  profiles  hare  been  calculated  for  free-stream 
Mach  numbers  of  0 ,  1,2,  8,  and  o.  For  each  Mach  number . 
velocity  profiles  are  presented  for  both  a  linear  and  a  0.7  G- 
power  variation  of  viscosity  with  absolute  tern pe rat urt .  The 
calculations  for  a  linear  variation  an  much  simp! tv  than  those 
for  a  O ,i ti- poirt /’  variation,  it  <s  s/norn  iiou  og  sihcinnj  tin 
constant  of  proportionality  in  the  linear  approximation  such 
that  it  gins  the  correct  calm  for  tht  viscosity  in  tht  high- 
temperature  part  of  the  mixing  layer ,  the  resulting  velocity 
profiles  are  in  excellent  agreement  with  those  calculated  by  a 
0.7  G- power  variation. 

INTRODUCTION 

The  velocity  profile  for  turbulent  mixing  at  constant 
pressure  of  an  incompressible  stream  with  a  dead-air  region 
has  been  calculated  by  several  investigators,  principally 
Tollmien  (reference  1).  These  calculations  agree  well  with 
the  available  experimental  data,  although  the  conventional 
assumptions  regarding  the  mixing  length  of  a  turbulent 
flow  have  since  been  shown  by  experiments  to  be  incorrect 
(reference  2).  The  many  difficulties  encountered  in  making 
precise  turbulent-mixing  calculations  are,  of  course,  a 
consequence  of  the  extremely  complicated  mechanism 
govering  all  turbulent  flows.  Jn  contradistinction  to  the 
case  of  turbulent  mixing,  the  mechanism  involved  in  laminar 
mixing  is  relatively  simple,  and  the  inalheinaiieal  relation 
between  stresses  and  velocity  gradients  for  laminar  flow 
is  well  kmwi>.  The  velocity  profiles  for  laminar  mixing, 
however,  apparently  have  not  as  yet  been  calculated  even 
for  the  case  of  incompressible  flow.  It  is  the  purpose  of  the 
present  paper  to  calculate  the  velocity  profiles  for  laminar 
mixing  (starting  with  zero  boundary-layer  thickness)  of 
an  air  stream  of  arbitrary  temperature  with  a  dead-air 
region  also  of  arbitrary  temperature.  In  eases  where  a 
laminar  boundary  layer  of  appreciable  thickness  exists  at 
the  point  where  mixing  begins,  the  results  given  herein  are 
not  directly  applicable  in  the  initial  part  of  the  mixing 
region.  For  such  cases,  it  is  necessary  to  make  some 
supplementary  approximation  in  order  to  apply  the  results. 

Since  the  practical  applications  of  laminar-mixing  phe¬ 
nomenon  usually  involve  the  flow  of  a  gas,  the  present 

S78IMS—  '><! 


1  analvsis  includes  tin*  effects  of  eomnressihilit  v  Kxainples 
ol  typical  flows  wherein  laminar  mixing  occurs  can  hr  lounu 
in  the  flow  of  small-scale  jet  pumps,  in  the  flow  behind  the 
intersection  of  shock  waves  of  unequal  strength,  and  in 
the  flow  immediately  behind  the  base  of  a  body  which  has 
a  laminar  boundary  layer.1 

SYMBOLS  AND  NOTATIONS 

('  constant  of  proportionality  between  viscosity  and 
temperature 

cv  specific  heat  at  constant  pressure 

k  coefficient  of  boat  conduction 

!.  ,  ham*  :  ic  !eng:!; 

M  Mach  number 

p  static  pressure 

Pr  Prandtl  number 

5  Sutherlands  constant,  approximately  216°  F  for  air 

T  absolute  temperature 

L7.  free-stream  velocity 

u,  v  velocity  components  in  x,  y  directions,  respectively 
x,  y  Cartesian  coordinates 

6  thickness  of  mixing  layer,  taken  between  points  where 

the  velocity  is  0.01  and  0.99  of  the  free-stream 
velocity 

a>  exponent  of  viscosity  variation  with  temperature 
p  mass  density 

y  coefficient  of  viscosity 

v  kinematic  coefficient  of  viscosity 

\f  stream  function 

f  dimensionless  independent  variable 

SUBSCRIPTS  AND  SUPERSCRIPTS 

*  dimensionless  variables  as  defined  in  equation  (14) 

®  free-stream  conditions 

o  stagnation  conditions  of  the  free  stream 
d  conditions  in  the  dead-air  region 

BASIC  EQUATIONS  AND  ASSUMPTIONS 

BASIC  EQUATIONS 

A  schematic  illustration  of  the  flow  under  consideration  is 
shown  in  figure  1.  In  order  to  make  the  laminar-mixing 
process  amenable  to  calculation,  the  usual  assumptions  are 
made  that  the  layer  affected  by  viscosity  is  thin  and  has  zero 

i  The  present  analysis  was  undertaken  as  part  of  an  investigation  of  this  latter  problem, 
and  originally  appeared  as  Appendix  B  of  a  thesis  "Base  Pressure  at  Supersonic  Velocities,” 
submitted  to  the  California  Institute  of  Technology.  June  1WK.  The  results  of  some  supple¬ 
mentary  computations  not  given  in  the  thesis  have  been  added  for  sake  of  completeness  in 
the  present  report. 
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pressure  gradient.  Under  these  conditions  the  formal  pro¬ 
cedure  for  estimating  the  order  of  magnit  ude  of  the  various 
terms  in  the  complete'  Na vier-Stokes  equations  for  viscous 
compressible  How  can  he  carried  through  in  precisely  the 
same  manner  as  is  done  in  the  classical  (Brandt  1)  treatment  of 
laminar  boundary-layer  flows.  In  so  doing,  the  dynamic 
equation  for  the  x  direction  reduces  to  the  familiar  boundary- 
layer  momentum  equation 


bv  .  bu 
pudx+prdy 


b  /  bu\ 
by  v  tyf) 


(1) 


while  the  dynamic  equation  for  the  y  direction  reduces  to 
zero  on  both  sides.  In  passing  from  the  Navier- Stokes  equn 
tion  to  tl)  it  is  to  be  noted  that  the  usual  boundary-laver 
assumption  6/j<<1  is  violated  in  th  eimmediate  vicinity 
of  point  0  (fig.  1)  just  as  in  the  case  of  boundary-layer  flow 
over  a  plate. 

By  employing  the  same  considerations  on  order  of  magni¬ 
tudes  as  were  used  for  the  complete  Navier-Stokes  equations, 
the  complete  differential  equation  representing  the  balance 
of  energy  in  viscous  compressible  flow  reduces  to 


d,.T  s  0  .  .  7 

in  TP'"  d,/ 


which  is,  of  course,  the  usual  energy  equation  for  laminar 
boundary-layer  flow.  In  addition  to  equations  (1)  and  (2), 
the  equation  expressing  conservation  of  mass  is  needed  : 


5  (aw)  ,  d(p0)_ 
bx  by'~ 


For  a  given  gas  the  variation  of  m  and  cv  with  temperature  is 
known;  hence,  the  foregoing  system  of  three  partial  differen¬ 
tial  equations  is  completed  by  the  addition  of  the  equation  of 
state  for  a  region  of  constant  pressure 


T 

T 

^  OD 


P® 

P 


(4) 


ASSUMPTIONS 


In  order  to  solve  the  above  system  of  equations,  the  follow¬ 
ing  assumptions  are  made: 

1.  rr  =  constant 

2.  rr~cPiijk—  1 


3. 


’  where  C  is  a  constant  depending  on  T 


and  Tm 


The  second  of  these  assumptions  if  often  made  in  calculating 
boundary-layer  flows  when  only  the  velocity  profile  is  desired 
and  not  tin*  thermal  characteristics.  The  difference  between 
the  boundary-layer  velocity  distributions  for  Pr~  1  and 
Pr= 0.73  is  small,  as  is  clearly  shown  by  the  numerical  results 
of  Emmons  and  Brainerd  (reference  3).  Since  the  mixing- 
layer  and  boundary-layer  flows  differ  only  in  the  boundary 
conditions  and  not  in  the  differential  equations,  the  effect  of 
assumption  2  may  be  expected  to  be  similar  in  the  two  types 
of  viscous  flow.  At  moderate  supersonic  Mach  numbers,  the 
use  of  Pr~  1  for  air  does  not  introduce  more  than  1-  or 


2-percent  error  in  the  boundary-layer  velocity  profile;  and 
hence,  for  all  practical  purposes,  the  mixing-layer  velocity 
distribution  calculations  for  Pr=  1  should  be  sufficiently 
accurate  for  air. 

Assumption  3  needs  some  explanation  since  the  introduc¬ 
tion  of  a  constant  C  differing  from  unity  in  the  approximate 
relation  between  viscosity  and  temperature  apparently  has 
not  been  used  in  previous  work.  Usually  C  is  taken  as  unity, 
and  in  such  cases  the  approximation  m/m=°  =  (T/Tm)w  gives 
reliable  results  for  a  fixed  provided  the  free-stream  tem¬ 
perature  is  restricted  to  a  certain  ramre.  By  introducing  the 
lactor  (  ,  the  approximating  equation  can  m*  inauc  to  give 
the  same  value  as  a  more  exact  equation  at  any  desired 


Figure  X.— Schematic  drawing  of  the  flow. 


temperature  in  the  mixing  layer  regardless  of  71*  or  w. 
Assuming  that  Sutherland’s  equation 


m  /TV'2  Ta+S 

Moo  V  T®/  T+S 


(5) 


represents  the  true  variation  of  viscosity  with  temperature, 
then  the  approximate  equation  can  be  made  exact  at  any 
given  temperature  T  by  means  of  the  relation 

In  particular,  if  the  approximating  equation  is  linear  in 
temperature  (w=  1 ,  thereby  greatly  simplifying  the  boundary- 
layer  equations)  and  the  viscosity  is  matched  at  the  tem¬ 
perature  I'd,  than  the  above  equation  gives 

r_  fflT.+S  m 


as  the  value  of  the  constant  (\  By  selecting  C  in  this  man¬ 
ner,  rather  than  taking  it  as  unity,  a  linear  variation  of 
viscosity  with  temperature  then  becomes  an  accurate  ap¬ 
proximation  in  the  inner  part  of  the  viscous  layer,  rather 
than  in  the  outer  part  where  the  viscous  stresses  are  less 
important . 

SOLUTION  TO  BASIC  EQUATIONS 

As  was  first  pointed  out  by  Prandtl  in  reference  4,  and 
later  used  to  advantage  by  Busemann  and  Crocco  (refer¬ 
ences  o  and  6,  respectively),  the  consequence  of  the  assump¬ 
tion  Pr—  1  when  applied  to  boundary-layer  flow  is  that  the 


temperature  becomes  a  fund  ion  only  of  t lit'  veloeit  v.  lienee 

cpT—j[n)  (8) 

By  substituting  i  bis  relation  into  equation  (2i  ami  using 
equation  (lj  in  conjunction  with  the  assumption  /V-  i,  it 
follows  that  the  energy  equation  is  automatically  satisfied 
if  t  he  function  f(u)  satisfies  the  ordinary  different ial  equation 


LAM  I  N  Ali  M  I  \  I  \ ( ;  ni  A  COM  I’K  KSS1 BLL  I  LI’ ID 
1 1  follows  t  hat 


di+ 1 ; 


=  0 


1  til  eg  I ‘;i  i  i  I  i  l  lil>  etjwi:  i  Jolt .  -  i. 

T—  T„  for  m  =  l  \ 
T-  T,}  for  u  —  0 


gives 


f  (  m"i  -----  r 


r  T~  <‘r  T,!  —  ”  +  "  [c,  ( T  —  T<1 }  ^  2  J 


p« 

~  p  dy} 


d\p 

dx 


(12) 


d^d*  d  .  ds  0  __  pr  0 

dx  dx  d\p^~  dx  ds  p«,  d\p 


On  0 
0/  ds 


Setting  s~x. 


ds  ds 

N  =0  and  -  =  1 
du  dx 


that  the  transformation  formulas  are  2 


0 

0/ 


pu  d 

Pec  d\f' 


(9) 


(10) 


(1 1) 


/  d  \  __  pr  0  ,  /  d  \ 

\dx),  poc  ”r  \  dx/, 


as  the  relationship  between  velocity  and  temperature. 
Since  the  temperature  determines  tin*  density,  equation  (11) 
also  provides  a  means  for  calculating  the  density  as  a  func¬ 
tion  of  the  velocity. 

Following  the  method  first  given  for  incompressible  flow 
by  von  Mises  (reference  7)  and  later  list'd  for  compressible 
How  by  von  Kantian  ami  Tsien  i reference  Si,  a  transforma¬ 
tion  is  made  to  a  new  set  of  independent  variables  (/.  i/0, 
whore  $  is  the  stream  function.  By  using  \p  as  one  of  the 
independent  variables,  the  continuity  equation  (3)  is  iden¬ 
tically  satisfied,  and  the  velocity  components  are  given  by 


Since  the  requirements  of  conservation  of  energy  are  ful¬ 
filled  by  equation  (11).  and  conservation  of  mass  by  (‘filia¬ 
tion  (12),  the  only  equation  now  remaining  to  be  satisfied 
is  the  momentum  equation  (1).  If  a  transformation  wore 
made  to  a  completely  new  set  of  independent  variables 
(.v,  \p).  thi*  transformation  formula  would  he 

d  _  d\p  d  ds  0  pu  0  ds  d 
dy~~ dy  dy  dy  ds  p*  d\p'  dy  ds 


and 


becomes 


(du 

\  .  / 

du  \ 

( 

du  \ 

\dx 

)  +P-'( 

dy  )r 

-  pu  [ 

dx) 

1  \ 

P" 

(  °U  \ 

it 

\  £>y), . 

a 

P« 

WA 

(*ntu 

m  efjiiation 

(l)  in 

the 

bn  d 

/ 

p  bv' 

P« 

dx  :  dt 

(  Up 

p(t  d\p  > 

) 

(13) 


This  can  be  put  in  dimensionless  form  by  introducing  the 
variables 

T 

Y„ 


p  — 


x 

L 

p 


7’*- 


=  CT*U 


svM~LC 


(14) 


Kxcept  for  the  parameter  C  appearing  in  the  definition  of  \p* 
and  p*,  these  variables  are  the  same  as  those  used  by  Karin  an 
and  Tsien  (reference  8).  Remembering  that  To*,  the 
free-stream  total-temperature  ratio,  is  given  by 


To 

t; 


7—1 

2 


tin'll  tin1  relation  (equation  (11))  between  temperature  and 
velocity  can  be  written  as 


T*  —  T(i  *  — 


7-1 


(15) 


The  momentum  equation  (13)  becomes,  using  p*T*=  1, 


dv 11 
dx* 


d 

dt 


(* 


*7’*“’  1 


(16) 


This  is  the  basic  equation  which  must  be  solved.  The 
boundary  conditions  of  the  problem  are  such  that  no  bound 
ary  layer  exists  at  the  point  whore  mixing  first  begins. 
l:ndor  these  conditions  the  velocity  profiles  will  be  similar 
at  all  points  downstream  of  the  origin,  and  hence  the  velocity 
u *  will  be  a  function  only  of  some  dimensionless  variable  f. 
This  dimensionless  variable  must  involve  both  \p*  and  /*, 
and  must  he  zero  at  the  origin  of  coordinates  since  the  mixing- 
layer  thickness  is  also  zero  at  the  origin.  Therefore,  let 

£~\f/*aX*h 

where  (i  and  b  are  pure  numbers  which  must  be  determined 
by  the  condition  that  both  sides  of  equation  (16)  for  v*  are 
functions  only  of  the  single  variable  f.  Setting 


:  Tin*  variables  held  constant  in  a  ditTerentiation  process  arc  cxpliciily  indicated  in  those 
cases  where  ambiguity  could  result  if  the  subscript  notation  were  not  used. 


</(f)  =u*T*u 


(17) 


sisms — r.n — 
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then  the  right  side  of  equation  (Hi)  can  be  written  as 

d  /gdu*  df\  ?>  (ndu*  ir,b\ 

dtp*  \  (is  <ty*/  0^"  \4  (I±  } 


K 


dt*V  d$ 


g^jr  a**a-') 


from  which  it  is  obvious  that  in  order  for  the  right  side  of 
equation  (16)  to  be  a  function  only  of  f,  it  is  necessary  that 
o=l.  With  u*=u*( f)  and  a—  1,  the  partial  differential 
equation  (16)  reduces  to  the  ordinary  differential  equation 


(Iu* 

9  df 


Consequently,  in  order  for  the  entire  equation  to  be  a  func¬ 
tion  only  of  f,  it  is  also  necessary  that  b=—  hence 


*  =  ^ 

4  *  yr*  \U„VwxC 


(18) 


and  the  ordinary  differential  equation  for  the  velocity  dis¬ 
tribution  now  reduces  to 


f du*_d  /  <7u*\ 

2  d{  ~dA°  WJ 


(19) 


Equation  (19)  is  the  same  differential  equation  that  was 
obtained  for  boundary-layer  flow  in  reference  8.  It  is  a  non¬ 
linear  differential  equation  since  g  depends  upon  the  velocit  y 
u*.  This  equation,  however,  can  easily  be  converted  into 
an  integral  equation  which  can  be  solved  by  the  method  of 
successive  approximations.  The  conversion  is  made  by 
temporarily  assuming  that  g  is  a  known  function  of  f  (instead 
of  u*)  and  formally  applying  the  standard  methods  for 
solving  first  order  linear  differential  equations.  The  result  is 

U*  =  (\  (Hrff+C,  (20) 

j  o  y 

where 

-  fr  f  «r 

F=e  J”  5  (21) 


The  boundary  conditions  are 

u*~- 1  at  ^  —  00 
v*  —  0  at  f  =  —  < 


(22a) 
1 22  b) 


Letting  uu*  be  the  value  of  v*  at  f  =  ().  equation  (20)  can  he 
written  as 


'.r  \j 


The  constant  (\  must  satisfy  two  requirements 


1  — w„ 

*00  Ji 

»/d  y 


( \  ----- 


dt 


•u  t 

g 


d  f 


( 2o » 


(24) 


Equation  (2;>)  is  an  integral  equation  for  ?/*.  since  both  / 
and  0  are  funet  ions  of  ?/  *.  By  simply  est  ima  t  imr  a  reasonable 


solution  for  u*  as  a  function  of  J,  a  first  approximation  to 
the  true  solution  can  be  calculated  from 

iw*=c,J»  i  ^■°uu* 

The  zero-order  approximations  0F  and  o9  can  be  calculated 
directly  from  ntz*  by  using  equations  (15),  (17).  and  (21). 
If  this  process  is  repeated  untill  a  given  approximation  is  the 
same  as  the  preceding  one  (to  the  degree  of  accuracy  desired), 
and  equation  (24)  is  simultaneously  satisfied,  then  the  solu¬ 
tion  to  ihn  oroUnm  L  obtain nA  The  itenitioii  vvoc<'«  turn< 
out  to  he  rapdiiy  convergent,  requiring  two  or  three  itera¬ 
tions  to  obtain  the  function  «*(£)  accurate  to  within  1  per¬ 
cent,  and  about  four  or  five  iterations  to  obtain  w*(f)  accu¬ 
rate  to  within  n  few  tenths  of  1  percent. 

In  order  to  change  the  function  u*(£)  back  to  the  physical 
coordinates  (j,  y)  a  simple  quadrature  is  necessary.  B\ 
definition  of  the  stream  function, 

(sw® 


dx+  \  U„v„xC d{=~  dy—-~  dx 

X  poo  P® 


Hence,  with  x  held  constant ,  integration  gives 


(25) 


from  which  ti/L7.  as  a  function  of 


V\  vx  cnn  ^^ermined. 


It  is  to  he  noted  that  no  graphical  or  numerical  differentia¬ 
tions  are  needed  at  any  point  in  the  above  iteration  process, 
only  quadratures  are  required. 

As  is  evident  from  equation  (16),  the  assumption  m *~CT *, 
that  is,  1,  makes  the  momentum  equation  (16)  independ¬ 
ent  of  tempeiature,  and  hence  density.  Consequently,  with 
o)~  1,  the  solution  to  equation  (16)  in  ( x ,\fr)  coordinates  is 
independent  of  Mach  number.  For  zero  Mach  number, 
T*=  1,  ('=  1 ,  and 


Lsmg  this  relation  the  solution  in  physical  (/,?/)  coordinates 
for  a  linear  variation  of  viscosity  with  temperature  is  obtained 
from  the  solution  for  zero  Mach  number  by  substituting 
equation  (15)  into  equation  (25).  This  gives 


!  I  00  rji  *  /  '  f  00  \ 

?/\  vZxC  Id  \  ,J  \  r„xL.„ 

‘  MJ  I  ‘  v*dt-r(Tb*-rS)i  (2«> 

Jo 

where  the  integral  in  the  second  term  on  tin*  right  side  is 
carried  out  for  //*(f)  corresponding  to  Mx  -  0. 


LAMINAR  MIXING  OK  A  (’OMPliKSSIHLK  KLT’ID 


RESULTS 

Numerical  calculations  of  the  velocity  distribution  have 
been  made  for  the  follo\viii«r 

1.  n*  =  T*°-n;  A/»  =  0,  1,  2,  3,  and  5 

2.  fi*  =  CT*;  M„  =  0,  1,  2,  3,  and  5 

The  various  solutions  for  case  2  are  obtained  directly 
(equation  (26))  from  the  solution  for  A/»=0  without  carry¬ 
ing  out  the  laborious  iteration  process  that  is  necessary  to 

•  l  for*  . . . . ,lt .  M  I, I 

coordinates  (y^  Um/y»x  as  independent  variable)  liave  been 
calculated  for  the  case  Ta—T{i.  If  the  dead-air  tempera¬ 
ture  is  radically  different  from  the  free-stream  stagnation 
temperature,  the  proper  velocity-distribution  curves  can  be 
obtained  by  carrying  out  the  integration  indicated  in  equa¬ 
tion  (2.V).  since  tin*  function  in  (xs\p)  coordinates  is 

independent  of  the  thermal  boundary  conditions  of  the 
problem. 

Curves  of  w*(f)  are  shown  in  figure  2  for  various  Mach 
numbers.  Tiie  corresponding  curves  in  the  physical  plane 
are  shown  in  figure  3  for  tin*  case  T*°-76,  and  in  figure  4 
for  the  case  v*  =  CT*.  In  the  latter  two  figures  the  familiar 
Blasius  curve  for  the  incompressible  laminar  houndary-layer 
flow  is  shown  for  purposes  of  comparison.  The  constant  C 
that  is  used  in  figure  4  is  determined  by  matching  the  viscos¬ 
ity  coefficient  at  the  temperature  Td*—Tv*  according  to 
equation  (7). 

The  particular  curves  shown  in  figure  4  apply  for  Tn  =400° 
R.  Curves  for  any  other  temperature  level  T„  differ  only 
in  the  constant  factor  C. 


CONCLUDING  REMARKS 

A  comparison  is  shown  in  figure  5  which  illustrates  the 
j£oo< I  agreement  bctwerui  vclocitv  distributions  calculated 
for  the  two  approximations,  a*  =  7’*° *70  and  n*  —  ('T *.  At 
a  Mach  number  of  2  or  less,  the  curves  for  a  linear  variation 
of  viscosity  with  temperature  virtually  coincide  with  the 
curves  for  a  0.76-power  variation.  For  general  use  the 
linear  approximation  is  recommended  since  it  gives  results 
which  are  practically  as  accurate  as  the  former,  yet  does  not 
require  a  laborious  iteration  solution  to  be  worked  out  for 
each  Mach  number 

in  general  the  laminar-mixing  layer  is  several  times 
thicker  than  the  laminar  boundary  layer,  as  is  illustrated  in 
figures  3  and  4  where,  for  purposes  of  comparison,  the  Blasius 
profile  is  also  shown.  The  rate  of  growth  of  mixing- layer 
thickness  with  increasing  Mach  number  is  somewhat  larger 
than  the*  corresponding  rate  of  growth  for  a  laminar  boundary 
layer.  The  curves  in  figure  if  indicate  a  value  of  roughly 

=  l+o.l  lil/.* 

o.v-0 

for  rate  of  growth  of  t ho  mixing  layer:  whereas  the  corres¬ 
ponding  value  for  a  laminar  boundary  layer  (see  reference 
9,  for  example)  is  approximately 

~  =l+0.0<Jil/„2 

6a/-o 

This  difference  is  to  be  expected  since  a  larger  percentage  of 
low-density  air  exists  in  a  mixing  layer  than  in  a  boundary 
layer. 


LAMINAR  MIXING  OF  A  COMPRESSIBLE  FLUID 


Figure  5.— Comparison  of  velocity  profile  for  a  linear  and  a  0.76-power  variation  of  viscosity 
with  temperature. 


The  foregoing  statements,  which  indicate  an  increase  in 
mixing-layer  thickness  with  increasing  Mach  number,  are 
based  on  the  assumption  that  the  Reynolds  number  (T’^rjv^ 
is  held  constant  while  the  Mach  number  is  varied.  In  most 
experimental  apparatus  the  Reynolds  number  changes  con¬ 
siderably  with  a  variation  in  Mach  number.  Consequently, 
depending  upon  the  particular  experimental  method  em¬ 
ployed,  the  observed  rate  of  mixing  in  the  z  direction  may 
be  either  increased  or  decreased  if  the  Mach  number  is 
increased. 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif.,  Jan.  f>,  191+9. 

REFERENCES 

1.  Tollmien,  Walter:  Calculation  of  Turbulent  Expansion  Processes. 

NACA  TM  1085,  1945. 

2.  Liepmann,  Hans  Wolfgang,  and  Laufer,  John:  Investigations  of 

Free  Turbulent  Mixing.  NACA  TN  1257,  1947. 

3.  Emmons,  H.  W.,  and  Brainerd,  J.  G.:  Temperature  Effects  in  a 

Laminar  Compressible-Fluid  Boundary  Laver  Along  a  Flat 
Plate.  Jour.  App.  Mach.,  vol.  8,  no.  3,  Sept.  1941,  pp.  A-105- 
A-110. 

4.  Prandtl,  Ludwig:  Eine  Beziehung  Zwischen  W&rme  austausch  und 

Stromungswiderstand  der  Flussigkeiten.  Phvsikalischc  Zeits- 
chrift,  vol.  11,  1910,  p.  1072. 

5.  Busemann,  Adolph:  Gas  Dynamics,  Theory  and  Concepts.  Chapter 

I,  Paragraph  9  in  Handbuch  der  Experimentalphysik,  vol.  4, 
Part  I.  Leipzig,  Akademische  Verlagsgesellschaft,  1931. 

6.  Crocco,  Luigi:  Transmission  of  Heat  From  a  Flat  Plate  to  a  Fluid 

Flowing  at  a  High  Velocity.  NACA  TM  690,  1932. 

7.  von  Mises,  Richard:  Bemerknung  zur  hydrodynamik,  Zeitschrift 

fur  angewandte  Mathematik  und  Mechanik,  vol.  7,  1927. 

8.  von  Kdrm&n,  Th.,  and  Tsien,  H.  S.:  Boundary  Layer  in  Compres¬ 

sible  Fluids.  Jour.  Aero.  Sci.,  vol.  5,  no.  6,  April  1938,  pp. 
227-232. 

9.  Lees,  L.:  The  Stability  of  the  Laminar  Boundary  Layer  in  a  Com¬ 

pressible  Fluid.  NACA  Rep.  876,  1947. 


U.  1.  GOVERNMENT  PRINTING  OFflCE:  l»10 


X 


Positive  directions  of  axes  and  angles  (forces  and  moments)  arc  shown  by  arrows 
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Rolling . 

L 
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V 
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Pitching _ ! 

Af 
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Pitch _ 

0 

V 

q 
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N 
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!  ip 

r 

Absolute  coefficients  of  moment 

r==.  L  rM 
1  f[bS  Vt  (jcS 

(rolling)  (pitching) 


qbS 

(yawing) 


Angle  of  set  of  control  surface  (relative  to  neutral 
position),  5.  (Indicate  surface  lyy  proper  subscript.) 


4.  PROPELLER  SYMBOLS 


D  Diameter 

p  Geometric  pitch 

plD  Pitch  ratio 

Inflow  velocity 
Si  i  p  s  t  re  am  v  e  I  o  c  it  y 

Thrust,  absolute  coefficient  CT 


1" 
T  T 
*  < 

T 

Q 


\ 


prrD* 

Torque,  absolute  coefficient  Ca— — 

prrD° 


c3 

v 

n 


Power,  absolute  coefficient  Cr>~ 


pn*D6 


Speed-power  coefficient  =  „ 

V  I  n* 

Efficiency 

devolutions  per  second,  rps 

Kifcctive  helix  angle = tan-1  ( ^ 
°  V2irr»/ 


1  hi)  — 70.04  kg-m/s-  ffi'Q  fi-lb.'sec 
1  metric  horsepower— 0.0803  hp 
1  rnph  =  0.4470  mps 
1  nip.?  — 2.2369  mpb 


5.  NUMERICAL  RELATIONS 


]  lb  =  0.4 53(5  kg 
1  kg=  2.2040  li) 

1  i n i --  POOO.Oo  m  =  5.2S0  ft 
1  m  =  :k2808  ft 


